An effective spin-orbital Hamiltonian for the double perovskite Sr2FeW06: 

Derivation of the phase diagram. 

S. Di Matteo,i'2 G. Jackeli,^ and N. B. Perkinsi'^ 

^ Laboratori Nazionali di Frascati, INFN, Casella Postale 13, 1-00044 Frascati, Italy 
'^INFM UdR Roma III, via della Vasca Navale, 00100 Roma, Italy 
'^Institut Laue Langevin, B. P. 156, F-38042, Grenoble, France 
^BLTP, JINR, Dubna, 141980, Russia 
(Dated: February 2, 2008) 

We formulate a superexchange theory of insulating double-perovskite compounds such as 
Sr2FeW06. An effective spin-orbital Hamiltonian is derived in the strong coupling limit of Hubbard 
model for d-electrons on Fe and W ions. The relevant degrees of freedom are the spins S = 2 and 
the three-fold orbital degeneracy of Fe^+-ions. W-sites are integrated out by means of a fourth- 
order perturbative expansion. The magnetically and orbitally ordered ground states of the effective 
Hamiltonian are discussed as a function of the model parameters. We show that for realistic values 
of such parameters the ground state is antiferromagnetic, as experimentally observed. The order 
found is of type-H, consisting of {111} ferromagnetic planes stacked antiferromagnetically. The 
orbital order energy scale found is one order of magnitude less than the spin one. 
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I. INTRODUCTION. 

Transition metal oxides have been the subject of nu- 
merous experimental and theoretical studies after the dis- 
covery of a variety of novel physical phenomena and a di- 
versity of ordered phases.^ Their ground state can be in- 
sulating, metallic or superconducting, depending on tem- 
perature and stoichiometry and their magnetic structure 
varies from ferromagnetic (FM) to different antiferromag- 
netic (AFM) types. The insulating behavior is gener- 
ally attributed to a large Coulomb interaction that local- 
izes the d-electrons in Mott-Hubbard or charge-transfer 
regimes.^ In this case low energy excitations can be de- 
scribed by the spin degrees of freedom through an effec- 
tive superexchange model. Often, the ground state of 
the transition metal ion possesses an orbital degeneracy 
in addition to its that originating from the spin. In such 
cases, orbital degrees of freedom are also incorporated in 
the superexchange theory and the system is described by 
means of an effective spin-orbital model. ^ 

The purpose of the present work is to formulate such a 
superexchange theory and derive a spin-orbital Hamilto- 
nian for the antiferromagnetic insulator Sr2FeW06. The 
experimental value of its Neel temperature Tn is quite 
low and reported in the literature as 16 — 37 Kii^ 
Sr2FeW06 belongs to a family of double perovskite com- 
pounds with general formula A2BB'06. It has, however, 
electronic and magnetic properties which differ signifi- 
cantly from those of the other members of the family. For 
example, Sr2FeMo06 and Sr2FeRe06, that have been the 
subject of the intense studies since the discovery of room 
temperature magnetoresistance in these compounds;^ are 
half-metallic ferrimagnets with fairly high magnetic tran- 
sition temperatures.^ One of the reason why Sr2FeW06 is 
so different from the other members of the group has been 
suggested by band structure analysis»i In this picture, the 
antibonding W(5d) states are pushed higher in energy by 



the stronger hybridization with oxygen p-orbitalsi^'^ (In 
the case of Re(5d) states this effect is compensated by the 
deeper 5d-level of Re ion.) Because of this, 5d-electron 
prefer to move away from W site, leaving it in a W^"*" 
{5(f) configuration, and stay on the Fe 3c?-level even at 
the cost of paying an extra Coulomb energy. As a re- 
sult, in Sr2FeW06 iron is in a Fe^^ (3d®) valence state^S 
and the moderately large Fe-W charge transfer gap, in 
the presence of strong Coulomb interaction on Fe sites, 
causes the insulating behavior. 

In the perovskite structure, the octahedral ligand field 
of oxygens splits d- levels into a low- lying triplet (t2g) and 
a higher doublet (e^). There are several experimental 
evidences^ suggesting that Fe^+ ion is in the high spin 
configuration tlg^g with S = 2. As four electrons occupy 
the threefold degenerate t2g-leve\s, the ground state of 
Fe^"*" has a threefold orbital degeneracy. This degeneracy, 
similarly to the spin one, can be described by pseudospin 
quantum numbers. The spin S = 2 and pseudospin t = 1 
are thus the building blocks of an effective theory for in- 
sulating Sr2FeW06. In the double perovskite structure 
of Sr2FeW06 Fe and W ions form a rock-salt structure. 
Each type of ions resides on a face-centered-cubic (fee) 
sublattice. The fee lattice is weakly frustrated and allows 
ordered anisotropic AFM structures. But weak frustra- 
tion of the lattice alone can not be responsible for the 
very low magnetic transition temperature. The reason 
why Tn is low, we believe, is as follows. The superex- 
change interaction, that couples magnetic moments of 
Fe ions, is mediated by two types of diamagnetic ions, 
along the path Fe-O-W-O-Fe, with 90° and 180° angles 
between Fe and W ions. Because of this, the resulting Fe- 
Fe exchange is rather weak and the corresponding transi- 
tion temperature is low. The nearest-neighbor (nn) and 
next-nearest-neighbor (nnn) Fe moments are coupled, re- 
spectively, by 90° and 180° Fe-W-Fe bonds. Depending 
on values and signs of the nn and nnn exchange energies, 
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different types of AFM structures can be stabilized on 
the fee lattice. To our knowledge no experimental de- 
termination of the kind of AFM order taking place in 
Sr2FeW06 has ever been undertaken. 

In this paper we derive an effective spin-orbital model 
from the Hubbard Hamiltonian of correlated Fe-3d elec- 
trons hybridized with uncorrelated 5d electrons of W and 
determine the magnetic ground-state structure as a func- 
tion of the model parameters. We assume that oxygen de- 
grees of freedom are integrated out and consider an effec- 
tive hybridization between Fe and W d-states. Perform- 
ing a fourth order perturbation expansion in the hopping 
parameters around the atomic limit we get an effective 
S = 2 Hamiltonian for a system with six electrons in 
five d orbitals. It may be worthwhile to note that our 
derivation is based only on the geometrical structure and 
electronic configuration of the system. As such it can 
be applied to a variety of systems presenting the same 
characteristics: any rock-salt with a high-spin f^gC^g con- 
figuration, like, e.g., FeO. 

The paper is organized as follows: In the next section 
we introduce the multi-band Hubbard Hamiltonian, per- 
form the fourth order perturbation expansion and write 
the effective spin-orbital Hamiltonian in which the or- 
bital degrees of freedom are described in terms of the 
pseudospin representation. The derivation of the effective 
Hamiltonian demands a lot of algebra with rather cum- 
bersome expressions. To facilitate the reader, we have 
tried to use, whenever possible, a pictorial representa- 
tion of states and processes, relegating both the fourth 
order perturbation scheme and all the final expressions in 
the Appendices A and B, correspondingly. In section III 
the classical ground state phase diagram of the effective 
spin-orbital model is discussed, as a function of model 
parameters, within a variational procedure. In the same 
section we also discuss the range of the external param- 
eters appropriate to the case of Sr2FeW06 and sketch 
some results. In section IV we draw our conclusions. 



II. EFFECTIVE HAMILTONIAN. 

Here we introduce the five-band Hubbard Hamiltonian, 
and perform a perturbative expansion in the hopping pa- 
rameters up to the fourth order to describe all the low- 
lying excited states of the Fe-W system in terms of spin 
and orbital degrees of freedom of iron sites, only. 

A. Ground state. 

We start from the multi-band Hubbard Hamiltonian: 
H = Hq + Hci + How + H[ , (1) 
where is the on-site term of the form: 

Ho = Ui'^nia]nial + -{U2- J) ^ riiaariia'a 
i.a i,(T,a^a' 



Ui and U2 are the Coulomb repulsion among electrons in 
the same and in different orbitals, respectively, and J is 
the Hund's coupling constant. Due to the cubic symme- 
try, the relation Ui = U2 + 2J holds. The annihilation 
and creation operators, diaa and d\^^ refer to Fe-3d or- 
bitals at site i, of type a (one of — y^, — r^, xy, 
xz or yz) and with spin a =t, [, n^aa = d]^^di„a- 

The octahedral crystal field splitting of d-orbitals into 
a Bg-doublet [x"^ — and — ) and a lower- lying t2g 
triplet (zy, xz and yz) is taken into account for both Fe 
and W ions by the termjiS 

i/c/ = E (riiaa +ni'aa) (3) 

where fii'aa counts the electrons on W sites i' . The sum- 
mation is for all i and i' running over Fe and W sublat- 
tices, respectively. Finally, the last term of the local part 
of the Hamiltonian, -ffoiv, represents the single particle 
energy level of the W ion with respect to Fe*^^^ : 

How = Act ^ n-i'aa , (4) 

z'qct 

where Act = E{Fe'^^^W^^^) - £;(Fe(6)w(°)) is the charge 
transfer energy, i.e., the minimum energy needed to re- 
move one electron from Fe^^-* ion and put it to an empty 
W'°^ site. Here and below an upper index denotes the 
number of d-electrons on a given ion. The kinetic term 
Hj. is considered as a perturbation: 

Hi = - E to^MaJr'^P + h.c) (5) 

The summation is over nearest-neighbor Fe-W sites, and 
all possible spin and orbital indices. Here d[i^^ {di'aa) 
creates (annihilates) a 5d-electron on W ions, with spin 
a in the orbital a. In a cubic lattice there is no hopping 
between eg and t2g. The transfer matrix tap is diago- 
nal in the t2g manifold for all the orbitals xy, xz, yz, 
with non-zero matrix elements t only for the bonds lay- 
ing in the corresponding plane (i.e., along x direction, 
the hopping integral is zero for yz orbitals and equal to 
t for the other two). Due to the shape of the Cg orbitals, 
their hybridization is different in the three cubic direc- 
tions thus leading to direction dependent hopping with 
the anisotropic transfer matrix elements tap {a = 'iz^—r^ 
and j3 — x"^ — y"^) given by: 

In constructing the effective Hamiltonian we also 
assume that the Hund's coupling J is strong enough 
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to win the competition with the crystal field splitting 
Agg. This implies that Fe^^-ion in the ground state 
is in the high-spin configuration ^1^6^ with total spin 
5 = 2, and justifies the exclusion from our zeroth 
order of the low-spin state, i^g^g {S = 0) and of the 
intermediate-spin state, i^g^g (-^ = 1); both higher in 
energy. We can write the ground state as a collection of 
atomic states \aj), each characterized by two quantum 
numbers: the spin S and pseudospin r. This latter 
describes the orbital occupation, namely the position of 
double occupied orbital in the t2g-manifold. Since there 
are three possibilities to put the double occupied orbital 
in the manifold, we can describe it by r = 1, with the 
same algebra of the usual spin operator. The orbital 
quantization axis can be selected arbitrarily, and we 
choose the following convention which pictorially can be 
presented as: 
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where \a) and \a') are two possible ground-states for Hq 
(all the Fe^'*'-ions in the ^Ig^g configuration, all 5d W- 
orbitals empty), while \l3m(i,n)) denotes all possible in- 
termediate excited states, as imposed by H^.. The first 
and third step states, |/3m(n))j rnust have five electrons 
on one Fe-ion and one electron in an adjacent W-ion 
(Pe(5)w^^)), while the second-step state, 1/3;), must be 
of the kind Fe^^^W^^^Fe^^) or Fe^^^W^^^Fe^^) on the three 
sites involved in the hopping process (see Fig. Our 
aim is to find out a representation of the effective Hamil- 
tonian, leading to this same eigenvalue equation, in terms 
of the spin and orbital degrees of freedom of iron only. 
We consider the full multiplet structure of the interme- 
diate states with five and seven electrons which can be 
grouped according to their energies and the total value of 
the spin. Moreover, depending on the orbital occupation 
of Fe(6) ion in the ground state and the orbital symmetry 
{eg or t2g) of the hopping electron, different intermedi- 
ate configurations can be reached. As already stated, in 
the cubic symmetry the hopping is allowed only within 
t2g-manifold and Cg-manifold, separately. 



In this notation, the upper boxes represent the Cg 
orbital doublet, Cg^'' = — (left upper box) and 
Cg^"* = — (right upper box). The lower boxes 
refer to the three-fold degenerate t2g orbitals. with 

tgg^ = xy (left box) , t'^J = yz (central box) and tfg = xz 

(right box). The atomic ground state energy is E^^ = 
I5C/2 — 8 J + 2 . Because of the three-fold orbital and 
the five-fold spin degeneracies, the ground state is 15- 
fold degenerate. Such a degeneracy is then lifted by the 
hopping term, as detailed below. 



B. Excited states and perturbation theory. 

In this subsection we briefly describe the perturbation 
scheme and classify the excited states according to their 
total spin S, their Sz projection and orbital configura- 
tions. Since we are interested in describing the superex- 
change interactions between iron sites, we need to go be- 
yond second order perturbation theory in the hopping 
Hamiltonian H[ and deal with fourth-order perturbative 
expansion. In fact, we neglect the direct Fe-Fe hopping 
and consider only the hopping through the W sites. The 
explicit form of H[ is such that the fourth order eigen- 
value equation simplifies to (see Appendix A): 



= (7) 

I 

Consider the excited states which involve t2g hopping, 

first. The intermediate atomic states with seven elec- 

(7'] 

trons FeJ have total spin S = 3/2. The manifold is 
12-fold degenerate (4(spin)x3(orbital)) and one of the 
possible states can be depicted as: 
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The three-fold orbital degeneracy arises from the t2g 

unpaired electron. The energy of this state is El = 
21^72 - 7J -f 2A. 

The intermediate states with five electrons split over a 
manifold with different values of total spin and energies. 

The excited states with maximum spin configuration 
have S — 5/2 and no orbital degeneracy. The 6- fold 
degeneracy is entirely due to the spin. The state with 
maximum Sz is: 



{a\H'^\p^){p^\H[\pi){l3i\H[\l3n){l3n\H[\a') 
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The energy is Ea = IOU2 — lOJ + 2A, the minimum 
excitation energy of Fe'^' states. 

All the other Fe'^^-' states with all orbitals singly occu- 
pied have total spin S — 3/2 and can be pictured as: 
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Each of these states is 4-fold degenerate because of the 
spin (only the states with maximum Sz are shown). They 
give the same contribution to the effective Hamiltonian 
and below, for convenience, we denote them with the 

(5) (5) 

common notation Fel . Their energy is i?^ — lOt/2 — 
5J-f-2A. 

Finally, there are the states with one unoccupied 
orbital. They have total spin S = 3/2 and the global 
degeneracy is 24 (4(spin)x6(orbital)). 
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The energies of these states are i?c+ = IOC/2 — 5 J- 



-2A 



and ErJ = IOC/2 -3J + 2A. 

Now we consider the intermediate states which can 
be reached by eg hopping. Those with seven electrons 
have spin 3 = 3/2 and can be presented as the following : 
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Ti T 



T Ti T 



where the normalization factors are defined as N± = 



These states have energy 



and A± = A±VJ2-H A2. 



(5) 

'd± 



IOC/2 - 4J + 2A T 
Vt/^ + A^ and their degeneracy is 4(spin) x3(position of 
doubly occupied orbital) x 2 (position of hole)=24. 



C. Classification of virtual processes and derivation 
of tiie spin-orbital Hamiltonian. 

In order to derive the effective spin-orbital Hamilto- 
nian in the simplest way it is useful to classify all the 
allowed hopping processes depending on the specific ge- 
ometry of the fee lattice. It becomes then possible to 
select the orbital dependence for each process and for 
a given projection of the spin coupling (whether FM or 
AFM) along the direction ij of the bond. The general 
form for the Hamiltonian turns out to be: 



(8) 



where ij denote the summation over both nn and nnn Fe- 



lons of the fee lattice. The terms O 



{F{A)) 



are the orbital, 



energy-dependent, contributions to the effective Hamil- 
tonian for FM or AFM bonds ij, respectively, which are 
explicitly given in Appendix B. The two spin-dependent 
terms are projectors over Si + Sj = A (FM coupling) and 
Si + Sj =0 (AFM coupling), respectively. Below we 
clarify the procedure adopted to derive the spin-orbital 
Hamiltonian (jSJ . We try to avoid as much as possible for- 
mulas and present the derivation in a pictorial way. Fol- 
lowing Fig. n we consider four types of processes (from A 
to D), according to the way the intermediate states can 
be reached. 

Depending on the kind and number of ions involved in 
the intermediate steps, it is possible to distinguish the 
following processes: 

A — processes : 
Fe(6)W(°)Fe(^) -> Fe^^^W^^^Fe^^) ^ Pe^^) W^'^^FeC^) 
Fe^^^W^i^Fe^s) Fe^^^W^o^Fe^^) 

B/C ~ processes : 
Fe(6)w(°)Fe(6) Fe^^^W^^^Fe^^' ^^(5)^^ (2)^^(5) 
Fe(^)w(i)Fe^^^ Fe^^^W^o^Fe^^^ 



Their degeneracy is 4(spin)x3(orbital t2g)x2(orbital 
eg)=24. Their energy is: eP = 2IU2 -7J + 3A. 

The states with five electrons arising from the eg 
hopping have also total spin 5 = 3/2 and can be drawn as 



D — processes : 
Fe(6)w(°)W(°)Fe(s) Fe^^^W^^^W^^^Fe^*^) 

Fe^^^W^^^W^^^Fe^'^) 



5 




Fe 



Fe 





FIG. 1: Different hopping processes: dashed and full lines cor- 
respond to different electrons. A- and B- processes exist for 
both 90° and 180° Fe-W-Fe angles (only 180° are shown). C- 
processes exist only for the 180°-angle and D-processes only 
for the 90°-angle. 



where the upper index denotes the number of electrons on 
each ion. A first basic distinction is easily drawn between 
the ^-process and all the others. The former involve the 
formation of a Fe^^-* state at the second step, while none 
of the latter type processes does. 

On the contrary, all B-, C- and D-processes are 
characterized by a Fe'^' state for both iron ions involved 
in a virtual process. The difference among them lies in 
the way these states are reached. For B- and C-processes 
there is only one intermediate W ion with two electrons 
on it, differently from the the I?-process, where two 
distinct W ions are involved. Finally, the difference 
between B- and C-processes lies in the behavior of the 
two virtual electrons. In the first case they go back 
to the original Fe-ion where they came from. On the 
contrary for C-processes there is an exchange of the 
two electrons, as shown in Fig. ^ and each of them 
moves toward the other Fe ion. For the sake of clarity it 
should be reminded that we are considering the specific 
case of the fee sublattice of Fe ions in Sr2FeW06: this 
allows some simplifications in the derivation of the 
effective spin-orbital Hamiltonian, -ffcff- In particular, it 
is possible to build the following table, where we list the 
contributions to Hcn depending on the characteristics of 
each single term. 
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AFM (Tl,nn) 



Each column in the table represents one of the four pos- 
sible processes drawn in Fig. ^ Each row represents the 
states of the Fe'^^-'-ion reached at the second step of the 



virtual process. For example, for B-, C- and D-processes, 
aa means that the virtual intermediate state is made of 

(5) (5) 

both FeJi , while ac indicates one FeJi -state and one 
Fef.^-state. In the case of the A process, only one Fe*^^^- 

state is present and "aa" refers to this single Fca -state, 
while ac has no meaning. The table gives us the magnetic 
character of the bond ij involved in the process (whether 
FM or AFM), if it requires a nn or a nnn contribution, 
and its orbital signature. In this respect Tl means that 
only one type of orbital is allowed in the process (ie, a 
contribution proportional to t^y), while T2 means that 
there must be two different orbitals (for example a con- 
tribution proportional to t'^yt'^^)- If there are no orbital 
or neighbors specifications, it means that both kinds are 
allowed. With this in mind, the derivation of the orbital 
part reported in Appendix B follows quite straightfor- 
ward. 



As an example, we show how to derive the term (jB4p 
(see Appendix B), corresponding to the A-process to- 
wards states c±. In order to have a state c± on a 
ion, the electron that jumps to the W-site and then forms 

(7) 

the FeJ state, must belong to one of the singly occupied 
t2g-states of the original Fe^^^^-ion. Only in this case, in 
fact, the Fe(^) -ion is left with one doubly occupied, one 
singly occupied and one vacant i2g-state. The electron is 
free to jump on W site, that is empty, and it has some 
restriction to hop to the next Fe^^^^ ion, because of the 
Pauli principle. Only some jumps are allowed and they 
do depend on the relative orbital occupations of the two 
Fe^^^-ions involved in the process. It is worthwhile to note 
that, if the spins of the two Fe(6) -ions involved were fer- 
romagnetically coupled, no hopping processes could have 
taken place because of the Pauli principle: the hopping 
electron would find all the singly occupied orbitals of 
the final Fe^^^-ion filled with the same spin projection. 
Thus, superexchange between c± states can take place 
only when the two Fe^®^ ions are coupled antifcrromag- 
netically, as shown in the previous table. 

Let us now analyze the orbital dependence. For sim- 
plicity we consider two specific configurations, of those 
that we called Tl and T2 in the table, e.g., the term 
depending on f^y and that depending on t^yt^^. The 
first hopping process is allowed only when t^^ ^ 1 and 
7^ 1, i.e. ccy-orbital states are singly occupied at both 
i and j Fe-sites involved. The second is active only when 
both and Tjz are either 1 or -1, i.e. either xy- or xz- 
orbital states are doubly occupied at both sites. During 
this process both orbital quantum numbers are either re- 
duced or increased by two, respectively. In the case of 
a doubly occupied xy-orbital, the xz-electron hops leav- 
ing a hole in this orbital state. The resulting state is 
not an eigenstatc of the Hamiltonian ^ and is then pro- 
jected over an eigenstate c± given by the superposition of 
states with xy- and a;z-holes. Therefore, an electron from 
the doubly occupied xy-orbital of Fe'-^-' can hop back to 
ion. This changes the initial orbital configuration 
Tiz = Tiz = 1 to the final Tiz = t^z = — 1 form. In terms 
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of the pseudospin quantum numbers, the two contribu- 
tions are: 
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where ± depends on which of the c± states are consid- 
ered. These terms enter into Eq. Q with a weight taking 
into account both spin (through Clebsch-Gordan coeffi- 
cients) and orbital projection of the single Slater deter- 
minant into the correlated state c±. The weight, or nor- 
malization coefficient (see the Table in Appendix B), of 
a given process is obtained by equating the matrix ele- 
ments of the Hamiltonian (jS)) to the corresponding ones 
obtained directly using the perturbation theory O. To 
calculate the weight of the proccess under consideration, 
for simplicity, we conider diagonal matrix elements of 
i?off in the spin subspace over the bond configuration 
SI = 2 and S] = -2. The SU[2) spin symmetry of the 
problem guarantees that consideration of the nondiago- 
nal matrix elements in the spin subspace, i. e. processes 
that involve spin flips, will lead to the same result. For 
the present process this weight is 7^^^ — ^ = ^{-^ ^ 



where the factor 1/8 normalizes \{Si ■ Sj 



the bond we consider and the remaining (-^ 



4)1 

\2 



8 for 



is due to 

the projection of the single Slater determinant into the 
states c±. The square power is due to the fact that the 
virtual state, not eigen-state of the local Hamiltonian, is 
projected twice into the corresponding eigenstate in Eq. 

O- 

All the other contributions, detailed in Appendix B, 
have been found in a similar way. 



III. GROUND STATE PHASE DIAGRAM. 

In this section we discuss the classical ground state 
phase diagram of the spin-orbital model, Eq. ©. We 
factorize spin and orbital degrees of freedom and per- 
form a mean-filed type analysis of the possible ordered 
phases. The factorization is justified because of the two 
different energy scales governing these degrees of free- 
dom: spins are coupled by virtual processes that involve 
both Cg and electron transfers, while only the latter 
electrons generate the coupling between the orbital de- 
grees of freedom, as Cg orbitals are half-filled. The ratio 
of spin and orbital exchange energies thus scales approx- 
imately as tg/f* and is much larger than unity in the 
situation of interest. Therefore, it is reasonable to aver- 
age out spin degrees of freedom, first, by considering all 
possible ordered magnetic structures on a fee centered 
cubic lattice and then for each type of magnetic order- 
ing find an orbital configuration, that minimizes the total 
energy by means of a variational procedure. 

We consider the following magnetic structures that are 
possible ordered states for a nn and a nnn Heisenberg 
Hamiltonian on a fee latticeii: 




FIG. 2: Type-II antiferromagnetic ordering of Fe moments on 
the face centered cubic lattice. 



i) Ferromagnetic state (FM); 

ii) Type-I antiferromagnetic (AFM-I) ordering, which 
consists of ferromagnetic {110} planes stacked antiferro- 
magnetically in the third direction; 

iii) Type-II antiferromagnetic (AFM-II) structure, 
consisting of ferromagnetic {111} planes, each being anti- 
ferromagnetically coupled to adjacent layers. This mag- 
netic structure is presented in Fig. 13 where only the Fe 
atoms are shown. 



A. Minimization procedure. 

The trial wave function can be written in general as 
follows: 



(9) 



where 1-0°) and \(j>f) refer to orbital and spin states at 
site i , respectively. Since H^g factors into an orbital 
H°g and a spin H^g part and its action involves two sites 
at a time, its average value over the above state takes the 
form: 



(10) 



The orbital part of the variational wave function can 
be written as 



^|0), + sin6',(cosVi|l)» + sini/;,| - 1)^) (11) 



allowing all three orbital states = 0)i, jr^ = 1)^ and 
\tz — — l)i with a relative weight determined by the min- 
imization procedure with respect to the variational pa- 
rameters 6i and ipi . The expectation values of H"^^ over 
the trial orbital wave function for both ferromagnetic and 
antiferromagnetic bonds are then expressed in terms of 
variational parameters 9i and and are quite lengthy. 
On the contrary, the average of spin variables is straight- 
forward in a mean field treatment. For a ferromagnetic 
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bond {Si ■ Sj + 6) = 10 and for an antiferromagnetic bond 
{Si ■ Sj — 4) = —8 in the classical ground state. 

The energy per Fe atom, E^^, is given by the sum of 
the energy contributions from all the bonds contained in 
the cell, divided by the number of atoms (each Fe-ion 
has 12 nn and 6 nnn). AFM-II phase has the biggest 
cell, with eight Fe atoms (see Fig.^J, and for generality 
we consider it as the unit cell for both FM and AFM-I 
phases, too. Hence, we have sixteen variational parame- 
ters (6*1 _8 and V'l-s for the eight sites of the unit cell). 
We have then performed a numerical minimization of the 
ground state energy with respect to all the sixteen inde- 
pendent angular variational angles 9i and ipi. By taking 
the absolute minimum we determine the orbital configu- 
ration that minimizes the total energy of the system for 
a given magnetic order. The ground state magnetic and 
orbital structure are determined by selecting the state 
with the lowest energy. The results are presented below. 



B. Model parameters 

In order to proceed further with the ground state phase 
diagram of HcS, we fix below the values of the parameters 
that have already been estimated in the literature and we 
try to give a reasonable guess for the remainings. 

• Crystal field splitting Ae^ : The value of the eg - 
^1.24 eV has been estimated for 
[l3 and we take this value in our 
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t2g splitting Ae^ 
Fe^"*" ion in Ref. 
analysis. 



• Exchange parameter J: The exchange part of 
Coulomb energy is known to be little screened 
in solids. It can thus be estimated trough the 
atomic values of Racah parameters B and G {J — 

+ C) . Given the values of B and C from Ref. 
Il3 . one obtains J ~ 0.78 eV. 

• Coulomb energy U2'- Intra-atomic Coulomb repul- 
sion is generally considerably screened in solids and 
can not be simply related to its atomic value. The 
estimates reported in the litera ture suggest that 
[/a = 5 - 7 eV (Refs. Il2ll3ll4j) . however we also 
consider variation of U2 in a wider range. 

• Transfer integrals t and t^,: We fix, in most of the 
cases, the ratio of Cg and t2g transfer integrals to 
be te/t = 3 as suggested by the Muffin-Tin Orbital 
theoryii^ We are not aware of any estimates of ab- 
solute values of the transfer integrals between Fe 
and W ions. However t2g hopping amplitude has 
been estimated to be < ~ 0.25 eV for Fe-Re bond.^^ 
We may consider the same value of t2g transfer for 
Fe-W bond too, since for both Re and W ions this 
transfer involves 5(i-shell. 

• Charge-transfer gap Act: We can estimate a lower 
boundary for the charge transfer gap between the 
high spin state of Fe'^^ and W ions in two possible 
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AFM-II 00-11 




FIG. 3: Ground state classical phase diagram of the spin- 
orbital model ||SJ in Act vs U2 parameter plane. Here AFM- 
II denotes type-II antiferromagnetic ordering, OO-I and OO- 
II correspond to two types of orbital ordering described in the 
text. 



ways as follows. First, the difference in Act for Re 
and W based systems can be approximately taken 
equal to the energy mismatch of 5o?-levels of these 
ions: AcT(W)-AcT(Re) = £;(5'i)(W)-£:(5d)(j^e)^ 
The latter mismatch is of the order of ~ 1.4 eVi^ 
A lower estimate of Act(W) is then obtained by 
putting Act (Re) = for the metalhc Re based 
compound. We get Act > 1.4 eV for the W based 
compound. Second, the insulating behavior of the 
system implies that the charge-transfer gap is larger 
than the bandwidth of those electrons that are 
involved in the lowest energy charge-transfer 
Taking the above reported value for t2g transfer one 
obtains the even more strict constraint Act > 2 
eV. 



C. Results and discussion 

In Fig. 13 we show the ground state phase diagram of 
the spin-orbital model ((SJl in Act — U2 plane. We fix the 
other parameters as discussed above and take te/t = 3. 
Because of this ratio, Cg superexchange, which is anti- 
ferromagnetic, overwhelms the superexchange due to the 
hg hopping. Therefore the ground state is always anti- 
ferromagnetic and it turns out to be of type-II. 

Given this magnetic structure, orbitals order in such a 
way to further reduce the ground state energy. The type 
of orbital ordering depends on whether 2Act ^ U2 (or- 
bital ordering of type-II, 00-11) or 2Act ^ U2 (orbital 
ordering of type-I, OO-I). In the former case virtual pro- 
cesses with singly occupied W-states (A-processes) be- 
come dominant. Instead in the latter case it is energeti- 
cally more favorable to create doubly occupied W states 
rather than to pay a Coulomb energy at Fe sites (B-, C-, 
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TABLE I: Types of orbital ordering with corresponding struc- 
tures 





OO-I 


OO-II 


OO-III 




> 


4- It^W /a/9 


\xz) 


2 


\xy) 


\xy) 


\xy) 


3 


\xz) 


\xz) 


\xz) 


4 


{\yz) + \xzy)/V2 


{\yz) + \xy))/V2 


\xy) 


5 


\yz) 


{\yz) + \xz))/V2 


\y^) 


6 


\xy) 


\xy) 


\yz) 


7 


\xz) 


\xz) 


\xz) 


8 


{\yz) + \xz))/V2 


{\yz) + \xz))/^/2 


\xy) 



D-processes). In Table |2 we list the types of orbital or- 
derings together with the corresponding structure. The 
numbers from 1 to 8 label the eight atoms of the unit 
cell shown in Fig. |21 The ket-vectors denote the type of 
doubly occupied t2g orbital at each atom for each orbital 
ordering. 

To understand the orbital structure of the AFM-II 
state it is convenient to divide the fee lattice of Fe ions, 
in four simple cubic sublattices. These four sublattices 
are generated by the pairs of nonequivalent atoms of the 
unit cell: (1,5), (2,6), (3,7), and (4,8) (see Fig.EJ. Each 
sublattice is AFM and the moments are coupled only by 
nnn exchange, while the coupling between sublattices is 
due to the nn exchange. If we denote by Ji and J2 the nn 
and nnn exchange constants of Fe spins, respectively, and 
consider only spin degrees of freedom, then the ground 
state energy in the AFM-II phase depends only on J2 , be- 
cause there is an equal number of nn ferro- and antiferro- 
bonds. Therefore the contributions from Ji to the ground 
state energy are canceled out. However, in the present 
case, due to the orbital part of the effective Hamiltonian 
inter-sublattice couplings (i.e. nn spin-orbital cou- 
plings) also contribute to the ground state energy. The 
antiferromagnetic superexchange for a given sublattice is 
maximized when t2g orbitals are ordered ferromagneti- 
cally and all equally populated. Switching on the inter- 
sublattice coupling does not affect the ferromagnetic-type 
orbital order within each sublattices. However, the type 
of occupied orbital state changes, in order to minimize 
the contribution of nn energy to the ground state. As 
a result, there are two types of orbitally ordered ground 
states (OO-I and OO-II), depending on the parameters 
(see Fig. 121. The structure of these orbital states is ex- 
plicitly given in Table 

We have also evaluated Ji and J2 averaged over the 
bond directions. The transition temperature can be 
directly linked to the nnn exchange energy J2 by the for- 
mula: Tjv — 2S{S + l)z J2/3 where z = 6 is a number of 
next-nearest neighbors. For S* = 2 and the experimental 
values of Tjy ~ 16 — 37 K one obtains J2 ~ 0.06 — 0.15 
meV. Fixing te/t — 3, U2 — 5 eV and varying Act in 
the range 2 — 7 eV we obtain the following range of J2: 
J2 = 0.03-0.14 meV for t = 0.20 eV and J2 = 0.08-0.35 




0.25 0.3 0.35 0.4 0.45 0.5 0.55 



te (eV) 



FIG. 4: Same as in Fig.|21in Act vs parameter plane. Here 
FM denotes ferromagnetic phase and OO-III corresponds to 
the orbital ordering described in the text. Filled and open 
circles stand for the FM-AFM-II phase boundaries for U2 = 5 
and 2.5 eV, respectively. 



meV for t — 0.25 eV. The exchange energy scales as the 
forth power of the transfer integrals and is thus very sen- 
sitive to the variation of t, nonetheless we can conclude 
that the superexchange theory gives the correct order of 
magnitude of the experimental transition temperature. 
It is worth mentioning that t2g electrons contribute ap- 
proximately only to one-tenth of the value of J2. Our 
analysis also suggests that Ji is antiferromagnetic as well, 
but around half of J2. Next-nearest neighbor spins are 
thus coupled stronger than nearest-neighbor ones. 

To complete the discussion on the ground state proper- 
ties of the spin-orbital Hamiltonian in the following 
we relax the constraint te/t = 3, fixing t = 0.25 eV and 
U2 — 5 eV, and consider the phase diagram in te and 
Act parameter plane. The result is presented in Fig. ^ 
For small values of eg hopping, a ferromagnetic phase is 
stabilized. Its orbital ordering (OO-III) is shown in Ta- 
ble In this configuration the system maximizes the nn 
ferromagnetic exchange, Ji > 0. The exchange between 
nnn neighbors is always antiferromagnetic, J2 < 0, and 
increases with ig. At some critical value of the Cg transfer, 
at which Ji + J2 — 0, the FM ordering becomes unsta- 
ble and a transition to the AFM-II state takes place (see 
filled circles in Fig. For smaller values of the Coulomb 
energy, U2 , the phase space for ferromagnetic ordering is 
enlarged and the transition to the AFM-II state occurs 
for larger values of t^- The FM-AFM-II phase boundary 
for J72 = 2.5 eV is shown as open circles in Fig. 0] 



IV. SUMMARY 

To summarize, we have developed a superexchange 
theory for insulating double-perovskite compounds such 
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as Sr2FeW06- The superexchange interaction, that cou- 
ples the magnetic moments of Fe^+ ions through the dia- 
magnetic ions, takes place on Fe-W-Fe bonds with 
90° and 180° angles. The theory has been implemented in 
an effective spin-orbital Hamiltonian formulated in terms 
of spin {S = 2) and orbital pseudo-spin (r = 1) degrees 
of freedom of the iron ion only, arising from the threefold 
degeneracy of i2c,-levels. We have studied magnetically 
and orbitally ordered states and constructed the classical 
ground state phase diagram of the effective Hamiltonian. 
We have shown that for realistic values of the model pa- 
rameters the ground state is antiferromagnctic, as experi- 
mentally observed. Our results suggests, that the antifer- 
romagnctic order is of type-II, i.e., {111} ferromagnetic 
planes stacked antiferromagnetically. This prediction can 
be tested by means of polarized neutron experiments. We 
have evaluated the exchange energies for these magnetic 
structure and showed that the theory gives an estimate of 
the transition temperature which is compatible with the 
experimental value. The magnetic structure is accompa- 
nied by an orbital ordering in which each cubic-sublattice 
of iron ions shows a ferro-type orbital ordering. The type 
of occupied orbital states changes from one sublattice to 
the other. The energy scale of orbital degrees of freedom 
is one order of magnitude smaller than the spin one. The 
is due to the fact that orbitals are coupled by t2g elec- 
trons, only, while both Cg (with larger transfer ampli- 
tude) and t2g virtual processes lift the spin degeneracy. 
Finally, we note that in the present paper we have just 
dealt with the electronic mechanisms to lift the ground 
state orbital degeneracy. 

However, there are two other possible mechanisms, not 



discussed here, that can lift the orbital degeneracy, such 
as the lattice distortion, or Jahn- Teller (JT) effect, and 
the spin-orbit coupling. The latter breaks the rotational 
spin invariance and opens a spin gap: Therefore a correct 
order of magnitude for this effect could be obtained by 
the experimental measurement of the spin- wave gap. As 
for the lattice mechanism, it should support the orbital 
ordering that minimizes the electronic energy of the sys- 
tem and thus enlarge the stability range of the orbital 
order, because the JT effect in this compound is mainly 
governed by the physics of isolated clusters rather than 
by cooperative effects. In fact, the oxygen octahedra sur- 
rounding Fe ions are independent one another, as they do 
not share any common oxygen-ion. Therefore JT distor- 
tions are not correlated as far as change in the electronic 
energy is concerned: Distortions of the octahedra around 
W ions with empty 5c?-level do not cause any considerable 
change in the electronic energy. 

In conclusion, if the electronic mechanism dominates 
the other two, then the orbital order can be stabilized 
only at very low temperatures, much lower than the mag- 
netic transition temperature and in between the system 
is magnetically ordered but orbitally disordered. In order 
to clarify which of the above mechanism is dominant fur- 
ther experimental studies of the magnetic and structural 
properties of the system are necessary. 
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APPENDIX A: DERIVATION OF THE FOURTH ORDER ENERGY CORRECTION. 



Consider the Hamiltonian: 

H = Ho + \H' (Al) 

where Hq is exactly solvable (HqiPq = Eo4'o) and A = is some "small" non-dimensional parameter. Following, for 
example R.efllfll it is straightforward to write the coupled equations to get eigenvalue and eigenvector corrections to 
£^0 and ipo up to the fourth order in A: 

(Ho - Eo)A - , (if' - E,)i^o + {Ho - Eo)^i = , (i/' - + [Ho - £o)V^2 - - , 

{H' - £;i)^2 + {Ho - So)^3 - E2A - E^A) = , 
{H' - Ei)^3 + (Ho - Eo)^Pi - E2^P2 - SaV^i - £^4^o = . 
Multiplying scalarly all equations by {ipo)* and solving with respect to the fourth order energy correction, E4, we get: 



P = ^fci-^im-^mn^nfc H^i H^^ S„^„ ^nfc H'kk / a 9^ 

' ;J^„^A;fe(^^-^')(^'^-^")(^'^-^") ^^{Ek-Ei){Ek-E^){Ek-E^) ^ ' 

H'klH'ikH'krnHlnk + H'kiHl^Hl^^^Hl,^ — H'l^iSlmHl^^H'i^f. 



{Ek - E,n){Ek ~ Ei)^ 

where we use notation: H'/.^^ = (?/'fc|i/'|i/'m)- Given the form of ifj, Eq. (Q, it follows that i/^^, — 0, ie, the second, 
fourth and fifth terms do not contribute. Also the third term can be neglected, because it does not give rise to an 
effective exchange between Fe-sites, being the product of two second-order processes. This justifies our formula in 
Sec.H.A. 
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APPENDIX B: ORBITAL TERMS AND ENERGY DENOMINATORS. 



Here we write explicitly the orbital contributions to the effective Hamiltonian ^ 
First we consider virtual processes wit 
be written in the following general form: 



First we consider virtual processes with t2s-electrons involved. All orbital contributions O^-J'^"^-'-' for a bond ij can 



0™'-E0":- where Oi^;p = ^Ol^;, 



IP nil 



p,w 



D'i 



(Bl) 



and Dp denotes an energy denominator, Ofj ^ is an orbital contribution expressed by means of the pseudospin 

operators, jp are normalization factors, which are given in the Table at the end of the Appendix B, P denotes one 
of the four virtual hopping processes A, B, C, and D, depicted in Fig.^ and I and I' are the kind of excited states 
involved in the hopping process (i.e. a, 5, and c introduced in Sec. II. B). 
In the following, we list all the terms, classified for the kind of process. 

j4-process: For processes involving excited states of a- and b- type the orbital contributions differ only because of 
the energy denominators: 



p.aa{bb) _ a 2 

— ^CT,a(f)) 



(7) , p(5) 



a{b) 



(B2) 



here and below Act,; denotes energy needed for an electron transfer from Fe*^^-' to W*^°^ when iron ion is left in a Fep-* 
state. This energy is given by Act,; — Act + ^ Ei^\ The state a is the lowest energy configuration of Fc^^^ and 
Act = Act, a is the lowest energy necessary for the Fe^^-* W*^"-* charge transfer. 
The corresponding pseudospin part is: 



O 



aa(bb) xy 



rl 



nMz + 1) 7^rj.{T,. - 1) + (1 - r/J + -^nAn. - l) -T,,(r,, + l) + (l - r,/,) + tUl - r^) 



-2^ 2 

3^ 



xy xz 



T- T- T- T- 



T- T- T ■ T ■ 

113 3 



tt.t 



xy^yz 



XZ yz 



(5) 

If the excited state is Fe^^ then the purely orbital contribution can be written as: 



(B3) 



/Sec _ .4 
^ij,A~ ''xy 



1 



(1 - rfz) + 7;nz{n, - 1) (1 - r/J + -r,,(r,, - 1) 



1 



(1 - 7-/J + T^nzinz + 1) (1 - r/J + -Tj,(t,, + 1) 



' yz iz jz 



± 



xy xz 



'l 'i 'j •] ^ 't 'i 'j •] 



± 



xy yz 



± 



'■xz'^yz 



(B4) 



where c can be either c+ or c_ , and respectively, either upper or lower sign should be taken in Eq. IB 41 Corresponding 
energy denominator is given by: 



^CT,c 



P(7) _^ £;(5) _ 2£;(6) 



(B5) 



B-process: It has only diagonal orbital contributions, as shown in the Table of section II. If aa-, bb-, ab- excited 
states are involved then the pure orbital part is the same, and can be written as: 



i3,B 



xy 



nzinz + 1)tjz{tjz + 1) + ^n,{nz ~ i)t,z{t,z - l) + i^,(i - t/J(1 - r/J 



(B6) 



For cc-hopping processes the orbital contribution Ofj^^ is equal to diagonal part of 0,f,^^ and reads as 



O 



+ t; 



(1 - rl) + ]^T,Mz - 1)] [(1 - rl) + \rj,Xr3z - 1) 
(1 - rl) + ir,,(r,, + 1)1 [(1 - t%) + \t,,{tj, + 1) 



' yz %z jz 



(B7) 
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As for ac- and be- processes, their orbital part are equal to the diagonal part of O^^^^^^ 



31 given by 



O 



ac{bc) 



xz / 
2 'iz 



The general form for B-processes energy denominators is 



(B8) 



D^^ = Act /A, 



CT,i^CT,i' 



A 



CT,i 



(B9) 



C-process: In this process two electrons with the same spin directions are exchanged between iron sites. The 
intermediate step involves the formation of a doubly occupied W site. This implies that two different orbitals must 
participate to the hopping process because of Pauli principle. The orbital contributions are thus only non-diagonal. 
They are the same for all processes of the kind aa, bb, cc, and ab: 



O 



II' 



I xy xz 



I I 3 3 



T- T- T- T- 

I t 3 J 



± 



tiJi 



'''i '''izTjzTj 



± 



(BIO) 



where the lower sign should be taken only when I = c_ and V 
energy denominators. 

For ac- and be- processes we have 



These processes differ one another only by the 



O 



ac(bc) 
i3,C 



I ''xy'^xz / + J- + -I- 

= ± [T T T T 

4 \ » 2 3 3 

I "xy^yz / + + _ _ N . 



xz^yz 



(Bll) 



where upper sign (+) and lower sign (— ) correspond to c = c+ and c = c_, respectively. 

We note that along 90° Fe-W-Fe bond (nn Fe ions) hopping is allowed only for one kind of orbital. Therefore, being 
orbitally nondiagonal, C-processes are effective only for 180° bond (nnn Fe ions), as shown in the Table of Sec. II. C. 
Energy denominators for all orbital contributions have the following common form: 



Z?" = Act, (Act, (' 



iCT,i 



iCT,i' 



(B12) 



D-process: This type of processes is effective only for 90° Fe-W-Fe bonds by definition [see Fig. ^ and thus can 
only be orbitally diagonal. All orbital contributions and energy denominators are equal to the corresponding ones of 
B-processes [see Eas. lB6IB9] . 

Now we consider Cg hopping processes. The contributions to the Hamiltonian from such processes does not depend 
on the orbital degree of freedom: 



dd 

(^dd _ IP 4.4 
'^ij,P - r)dd^e 
Up 



(B13) 



Here the extra factors, which depend on the type of process, and appear due to the anisotropic form of Cg transfer 
(0, are included in normalization factors given in the Table below. 

The energy denominators for different processes can be written in the following forms. For A-process 



r^dd 



^CT,d 



i5(7)+^(5)_2£;(6) 



(B14) 



and for B-, C- and D-processes 



D^Q — AcT,;AcT,i' Act,; + A 



iCT,i' 



(B15) 



where d and d! can be either or d^ excited states. 
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In the following Table we list all normalization coefBcients 7" : 





aa 


bb 


cc 


dd 


ac 


be 


ab 


A 


1/10 


1/10 


1/16 


l/8iV± 








B 


1/8 


2/25 


1/32 


l/(8iV|) 


1/20 


1/20 


2/25 


C 


1/10 


8/125 


1/40 




1/16 


1/25 


1/10 


D 


1/10 


8/25 


1/40 


1/(107V|) 


1/16 


1/25 


1/10 
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